.  AD-A149  048  ELIMINATION  OF  SENSOR  ARTIFACTS  FROM  INFRARED  DATA(U>  1/1 
NAVAL  RESEARCH  LAB  WASHINGTON  DC  R  STEINBERG  11  DEC  84 
NRL-8861 

UNCLASSIFIED  F7G  17/5  NL 


-  0 

END  1 

8m  aU  ■'■■'-  ‘  ^  ■«  1  *  ■"»*■■  ■*  1  **■  »*■»* 


microcopy  resolution  test  chart 

NATIONAL  BUREAU  Of  STANOAROS-I963-A 


Rf  PKODUCFfVAT  CiOVFRNMFNT’FXREflSE 


Elimination  of  Sensor 
Artifacts  from  Infrared  Data 

Richard  Steinberg 

Advanced  Concepts  Branch 
Optical  Sciences  Division 


NRL  Report  8861 


December  12,  1984 


Approved  for  public  release;  distribution  unlimiied. 


NAVAL  RESEARCH  LABORATORY 
Washington,  O.C. 


JAHl, 


85  01  08  158 


SEClm  CLASSIFICATION  Of  THIS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED _ 

2a  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 

I  lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 
4  PERFORMING  ORGANIZATION  REPORT  NUNISER(S) 


3  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution  unlimited. 

S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


NRL  Report  8861  _ 

6a  NAME  OF  PERFORMING  ORGANIZATION  1 6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

Naval  Research  Laboratory 

6c  ADDRESS  (City.  State,  and  ZIP  Code)  7b.  ADDRESS  (City.  State,  and  ZIP  Coda) 


Washington,  DC  20375-5000 

8a  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


8c  ADDRESS  (City,  Stare,  and  ZIP  Coda) 


1 1  TITLE  (Include  Security  Classification) 

Elimination  of  Sensor  Artifacts  from  Infrared  Data 

12  PERSONAL  AUTHOR(S) 

Steinberg,  R. 

13a  TYPE  OF  REPORT  |l]b  TIME  COVERED 

Interim  I  from _ ,  to  ___ 

16  SUPPLEMENTARY  NOTATION 


10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  I PROJECT  |T 

ELEMENT  NO  INO  It 


PROJECT 

TASK 

NO 

NO 

WORK  UNIT 
ACCESSION  NO 

EX38 0-280 


14  r>"E  OF  REPORT  (Year,  Month.  Day)  f  5  PAGE  COUNT 

1984  December  1 1  _ I  22 _ 


17 

COSATI  CODES  | 

FIELD 

GROUP 

SUBGROUP 

18  SUBJECT  TERMS  (Continue  on  reverie  it  necessary  , 

Infrared  BMAP 

IRST  Data  correction 


l  identify  by  block  number) 


19  ABSTRACT  (Continue  on  reverie  if  necessary  and  identify  by  block  number)  f 

v  r.'.— 

. —‘y  The  objective  of  the  Navy  Background  Measurements  and  Analysis  Program  (BMAP)  is  to  satisfy 

the  data  requirements  for  design  of  surface-  and  air-based  infrared  (IR)  search  and  tracK  devices.  The  IR 
backgrounds  data  are  acquired  with  a  sensor  that  has  high  spatial  resolution,  high  sensitivity,  and  simul¬ 
taneous  4  to  5dTm/8  to  lTunT capability.  Laboratory  tests  show  that  the  8  to  1  14tm,detector  array  has 
an  undesired  ground  resistance  that  introduces  weak  electrical  resistive  coupling  among  the  detector  chan¬ 
nels,  giving  rise  to  spatial  defects  in  the  long-wave  data.  Techniques  are  developed  for  inverting  the  focal 
plane  resistive'  network,  effectively  decoupling  the  detector  channels  and  eliminating  the  spatial  artifacts,  . 
without  smoothing  the  data  or  degrading  noise  performance.  £)ri«  -onW  -  v K>wah*  ( 

1  N  y  / 

y -  — '  -  .  1‘a^o.  I'dM^ct  ii)  n. 

fiv'r  rfrri/s  ^  l>  rtucrc£j  _ ^ _ 

20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

Bl/NqASSIFIEO/UNUMITED  □  SAME  AS  RPT  □  OTIC  USERS 

22*  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code?  22c  OFFICE  SYMBOL 

Richard  Steinberg  1  (2021  767-3045  Code  6520 

OD  FORM  1 473, 84  MAR  83  APR  edition  may  be  used  until  exhausted 


S3  APR  edition  mey  be  used  until  exhausted 
All  other  editions  ere  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


-  h 

nr 


CONTENTS 


BACKGROUND 


ARTIFACTS  APPEAR  IN  BAR-TARGET  SCANS .  2 


SUMMARY  OF  RESULTS 


FOCAL  PLANE  EQUIVALENT  CIRCUIT .  10 

BAR-TARGET  EXCITATION  .  12 

ITERATIVE  SOLUTION .  16 

CONCLUSIONS .  17 

ACKNOWLEDGMENTS .  17 


REFERENCES 


APPENDIX— Derivation  of  Equation  (35) .  18 


Aooaaalen  for _ 

"nTIS  ORAM  jf 

DTIC  TAB  O 

Unannounced  Q 

Justification—— 

By  . .... 

Distribution/ _ _ 

Availability  Codas 
™~"  Avail  and/or 
Dlst  Special 


ELIMINATION  OF  SENSOR  ARTIFACTS  FROM  INFRARED  DATA 


BACKGROUND 

The  objective  of  the  Navy  Background  Measurements  and  Analysts  Program  (BMAP)  is  to  satisfy 
the  data  requirements  for  design  of  surface*  and  air-based  infrared  (IR)  search  and  track  devices  [1,2]. 
The  BMAP  product  is  a  test  set  of  IR  background  images  for  use  in  off-line  simulation  of  alternative 
signal  processing  techniques  for  false-alarm  suppression.  The  sensor  used  to  obtain  the  background 
data  is  an  IR  scanner  constructed,  owned,  and  operated  by  Raytheon  Missile  Systems  Division,  Bed¬ 
ford,  MA.  Table  1  summarizes  the  measurement  system  characteristics. 


Table  1  —  Technical  Characteristics  of  the 
Raytheon  Dual-band  Scanner 


Pixel  size  (mrad,  square) 

1/3 

Elevation  channels/waveband 

16 

Azimuth  fieid-of-view 

2.2* 

Frame  rate  (frames/ s) 

1 

Word  depth  (bits/sample) 

12 

Sample  factor  (samples/dwell) 
NEI  (w/cm1,  array  average) 

3.44 

3.0-4.8  fim 

2.0  x  10"M 

7.6-11.3  pm 

1.0  x  iru 

The  scanner’s  long-wave  array  consists  of  two  columns  of  detectors,  separated  by  a  gold  common 
deposited  on  the  HgCdTe  detector  material.  Ground  connections  are  brought  out  from  the  gold  com¬ 
mon  both  from  the  top  and  from  the  bottom  (Fig.  1).  The  long-wave  array  has  IS  optically  active  ele¬ 
ments,  the  first  detector  channel  being  inactive. 


Fie<  1  —  Geometry  of  loot-wav*  photoconduc- 
tive  HgCdTe  detector  array.  Each  detector  fa 
square.  2  *  2  mil.  The  column  of  odd- 
numbered  detector*  i*  (operated  hom  the 
column  of  even-numbered  detector*  by  a  sold 
common  2  mil*  wide.  Ground  connection*  are 
made  at  point*  A  and  B  on  the  gold  common. 
Detector  No.  1  i*  optically  inactive. 
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The  detector  arrays  are  scanned  in  azimuth  over  approximately  a  2.8s  field  of  view.  As  each 
detector  traverses  its  2.8"  azimuth  sweep,  the  output  voltage  of  its  corresponding  postamplifier  is  sam¬ 
pled  S12  times.  Thus  each  long-wave  frame  consists  of  a  15  x  512  matrix  of  numbers,  where  the  small 
dimension  corresponds  to  detector  channel  or  elevation,  and  the  long  dimension  corresponds  to  time  or 
azimuth.  Of  the  total  2.8s  azimuth  scan,  about  0.6s  (the  first  100  samples)  is  filled  by  an  internal  cali¬ 
bration  pulse.  Thus  the  viewed  scene  external  to  the  sensor  is  actually  2.2s  in  azimuth.  We  omit  the 
calibration  pulse  in  displaying  the  data.  Censoring  the  first  112  samples,  each  long-wave  frame  then 
consists  of  a  IS  x  400  matrix  of  numbers. 


mM 


ARTIFACTS  APPEAR  IN  BAR-TARGET  SCANS 


Laboratory  bar-target  measurements  were  performed  in  the  interest  of  characterizing  the  transient 
response  of  the  Raytheon  dual-band  scanner.  Figure  2  shows  a  single  frame  of  long-wave  bar  target 
data  in  three-dimensional  (3-D)  perspective  format.  The  bar  was  oriented  perpendicular  to  the  direction 
of  scan  and  had  a  31.4°C  contrast  against  its  background.  The  data  displayed  in  Fig.  2  have  been  pro¬ 
cessed  by  NRL  as  follows. 


Sample  errors  introduced  by  the  digital  data  recorder  were  removed  [1]. 


•  The  data  for  even-numbered  channels  were  delayed  seven  samples  relative  to  the  data  for 
odd-numbered  channels  to  compensate  the  offset  between  even-  and  odd-channel  detector 
columns  (Fig.  1).  The  seven-sample  offset  follows  from  the  sample  rate,  3.44  samples 
per  dwell,  and  the  fact  that  the  gold  common  has  the  same  2-mil  width  as  the  detector 
elements  (Fig.  1).  The  distinction  between  the  true  two-dwell  offset  of  6.88  samples  and 
the  applied  integral  offset  of  7  samples  is  not  significant. 

•  Separate  offset  and  gain  parameters  were  calculated  and  applied  to  each  detector  channel 
to  compensate  detector  responsivity  nonuniformity. 

Before  inspecting  the  bar  target  measurements,  it  was  expected  that  the  preceding  sequence  of 
operations  would  constitute  the  major  part  of  data  reduction,  the  sole  remaining  step  being  a  single 
additional  gain  and  offset  correction  applied  to  the  entire  frame  for  radiometric  calibration.  This  expec¬ 
tation  was  overly  optimistic,  however,  as  seen  by  inspection  of  Fig.  2.  The  horn-shaped  spatial  artifacts 
in  the  data  are  seen  more  clearly  in  Figs.  3  and  4,  depicting  channels  2  and  3  of  the  data  frame  shown 
in  its  entirety  in  Fig.  2.  Data  for  all  even-numbered  channels  are  similar  in  appearance  to  those  of  Fig. 
3,  while  data  for  all  odd-numbered  channels  are  similar  in  appearance  to  those  of  Fig.  4. 

Raytheon  engineers  hypothesized  that  the  long-wave  data  artifacts  might  be  due  to  the  resistance 
of  the  ground  connections  to  the  HgCdTe  focal  plane  (Fig.  1),  introducing  a  small  amount  of  electrical 
resistive  coupling  among  the  detector  channels.  The  present  work  confirms  the  resistive  coupling 
hypothesis. 

SUMMARY  OF  RESULTS 

We  show  in  our  analysis  that  the  resistive  coupling  artifact  can  be  corrected  by  applying  the  fol¬ 
lowing  equation  to  the  long-wave  data: 


16 


En(t)  -  A~l  /„(*)  +  Bn  T  /.(f)  .  n  -  2,  3,  ...  16, 


(1) 


where  n  indexes  the  detector  channel,  t  is  time,  and  coefficients  A„  and  £„,  together,  are  30  correction 
constants.  The  15  waveforms  /„(/),  n  -  2,  3,  ...  16,  entered  as  input  to  Eq.  (1)  are  the  digitized 
postamplifier  outputs  in  need  of  correction  (e.g.,  Figs.  3  and  4  show  I2(t)  and  I3(t)  for  a  given  frame 


m3 


mm 
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Fig.  4  —  PosUmpUfier  output  waveform  for  LWIR  channel  3,  for  data  frame  ahown  in  its  en¬ 
tirety  aa  in  Fig.  J.  AH  odd-numbered  detector  waveforma,  3,  S, ...  15,  are  aimilar  in  appear¬ 
ance. 


of  data).  The  15  waveforms  £„(/),  n  -  2,  3,  ...  16,  are  the  corrected  data.  Equation  (1)  is  "instan¬ 
taneous,"  i.e.,  the  corrected  value  E„(tt)  is  generated  from  the  IS  instantaneous  samples  /j(ti),  /3(f|), 

. . .  and  does  not  depend  on  values  of  7„(f)  for  times  t  tj.  The  instantaneous,  or  "memory- 

less,”  property  of  Eq.  (1)  is  characteristic  of  resistive  coupling  as  opposed  to,  e.g.,  inductive  or  capaci¬ 
tive  coupling.  Since  Eq.  (1)  is  the  same  for  all  values  of  time,  we  generally  unburden  our  notation  by 
not  showing  explicit  time  dependence,  i.e.,  by  writing  the  waveform  quantities  E„(t)  and  I„(t)  as  E„ 
and  /„,  respectively. 

The  30  coefficients  in  Eq.  (1),  IS  values  each  of  An  and  B„,  are  determined  from  an  analysis  that 
requires  as  input  the  thirty  parameters  x„,  yn  obtained  from  one  frame  of  bar  target  data.  The  parame¬ 
ters  xm  and  yn  are  directly  measurable  from  the  original  waveforms,  as  shown  in  Figs.  3  and  4. 
Waveform  parameters  x„  and  y„  and  the  derived  correction  constants  A,1  and  B„  are  given  in  Table  2. 

The  same  bar-target  frame  in  which  the  resistive  coupling  defect  was  originally  manifest.  Fig.  2, 
serves  as  the  source  of  diagnostic  data  from  which  the  correction  constants  are  derived. 

The  effectiveness  of  the  correction  formula,  Eq.  (1),  is  illustrated  by  Figs.  S  to  7,  depicting  the 
corrected  versions  of  Figs.  2  to  4.  On  average,  the  resistive  coupling  artifact  is  reduced  tenfold  (Table 
3).  Apparently,  the  correction  is  excellent  when  applied  to  data  from  which  the  correction  constants 
are  derived.  An  obvious  question  is  whether  the  correction  is  data-  or  time-dependent,  i.e.,  whether 
the  correction  constants  derived  from  Fig.  2  will  serve  to  correct  data  having  different  varieties  of  struc¬ 
ture  than  Fig.  2,  obtained  some  time  subsequent  to  Fig.  2. 
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Table  2  —  Waveform  Parameters  and  Correction  Constants  for 
Long- Wave  Detector  Array.  Waveform  parameters  are  measured 
from  a  single  frame  of  bar  target  data  (Figs.  2 ,  3,  and  4).  Correc¬ 
tion  constants  are  used  with  Eq.  (1)  to  perform  data  correction. 


Detector  Channel  iWaveform  Parameter 


|  Correction  Constants 

4T* 

B„(x  irJ) 

1.01947 

0.68322 

1.01790 

0.63862 

0.99304 

0.99405 

1 

0.96181 

0.95269 

0.94931 

0.97060 

1.04890 

0.99043 

0.95794 

0.99591 

1.00078 

0.98810 

0.89753 

1.00925 

0.78908 

1.01933 

0.86301 

0.98741 

0.82757 

1.04250 

0.69041 

1.07017 

0.51001 

VV. 


RICHARD  STEINBERG 


E,(t) 


325  375 

SAMPLE  NUMBER 

Fig.  6  —  Channel  2  of  corrected  data,  extracted  from  Fig.  S 


E,0) 


325  375 

SAMPLE  NUMBER 

Fig.  7  —  Channel  3  of  corrected  data,  extracted  from  Fig.  S 
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Table  3  —  Performance  Summary  for  Resistive  Coupling 
Artifact  Correction  (1-bar  data).  Each  channel  waveform 
(e.g.  Figs.  3  and  4)  manifests  two  artifacts:  a  left  artifact 
and  a  right  artifact  The  relative  amplitude  of  each  artifact, 
left  and  right,  is  tabulated  as  a  "Percent  Error,"  calculated  as 
(100  yjx,).  Columns  labeled  "Avg"  are  the  numerical 
average  of  columns  labeled  "Left"  and  "Right*  The  column 
labeled  "Left"  under  "Original  Data"  can  be  calculated  from 
the  values  of  x,  and  ym  in  Table  2.  The  bottom  row  is  an 
average  over  channel  number,  i.e.,  an  array  average.  Thus, 
the  array-average  artifact  relative  amplitude  was  5.6%  in  the 
original  data,  but  is  reduced  tenfold  to  just  0.55%  in  the 
corrected  data. 


Detector 

Channel, 

a 


Percent  Error  (100  yjx,) 


Original 


Corrected 


Left  I  Right  I  Avg  8  Left  Right 


4.6  4.3  4.4  0.77  0.35 

4.8  4.2  4.5  0.50  0.06 

6.1  5.9  6.0  1.01  0.54 

5.6  5.6  5.6  0.57  0.67 

5.7  5.5  5.6  0.95  0.56 

6.8  6.3  6.6  0.70  0.25 

6.8  6.4  6.6  1.12  0.40 

7.2  6.7  6.9  0.74  0.33 

6.7  6.1  6.4  0.38  0.26 

6.7  6.2  6.4  0.69  0.23 

5.3  5.1  5.2  0.88  0.56 

6.6  6.5  6.6  0.68  0.64 

5.4  4.7  5.1  0.31  0.05 

5.4  5.1  5.2  1.24  0.35 


3.5 


Average  I  5.8 


3.3  fl  0.60  0.05 


0.74  0.35 


Figure  8  depicts  a  typical  even  channel  of  data  obtained  against  a  laboratory  target  consisting  of 
six  vertical  bars  of  varying  widths.  The  data  in  Fig.  8  were  obtained  in  measurements  made  about  three 
months  subsequent  to  the  Fig.  2  measurements.  The  correction  demonstrated  for  the  six-bar  data 
appears  almost  as  good  as  for  the  earlier  data.  On  average,  the  resistive  coupling  artifact  evidenced  by 
the  last  bar  of  the  six-bar  target  is  reduced  eightfold  (Table  4). 

We  note  that  the  later  data  were  corrected  with  correction  constants  derived  from  the  earlier  data. 

Figures  9a  and  9b  are  provided  as  a  final  example  of  how  the  resistive  coupling  artifact  appears  in 
original  and  corrected  field  data.  Figure  9a  shows  channel  2  of  a  long-wave  scan  across  the  moon, 
obtained  at  8:48  PM  on  14  August  1983  during  initial  sensor  field  trials  at  Montauk  Point,  LI.  The 
local  minimum  in  the  neighborhood  of  sample  #330,  which  is  due  to  the  resistive  coupling  effect,  is 
largely  eliminated  by  application  of  Eq.  (1)  to  the  data  (Fig.  9b). 
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Fig.  g  —  Channel  2  of  6-bar  data,  before  (/j)  and  after  (£})  resis¬ 
tive  coupling  correction.  Correction  constants  derived  from  1-bar 
data.  Fig.  2,  obtained  3  months  prior  to  the  6-bar  data. 


Table  4  —  Performance  Summary  for  Resistive  Coupling 
Artifact  Correction  (last  bar  of  6-bar  data).  Interpretation 
analogous  to  Table  3.  According  to  the  bottom  row  in  this 
table,  the  array-average  artifact  relative  amplitude  was  5.6% 
in  the  original  data,  but  is  reduced  eightfold  to  0.71%  in  the 
corrected  data. 


Detector 

Channel, 

n 

Percent  Error  (100  yjx„) 

Original  Data 

Corrected  Data 

Left 

Right 

Avg 

Left 

Right 

4.9 

mm 

4.3 

KY  TTCj 

■.933 

5.3 

w&m 

4.9 

1.13 

6.4 

5.1 

5.7 

1.19 

EES 

6.6 

6.1 

1.63 

0.64 

1.14 

mm 

5.4 

mi 

0.46 

0.70 

m 

6.5 

tfifl 

0.31 

0.54 

7.1 

6.5 

1.28 

0.46 

0.87 

9 

7.3 

WEm 

1.06 

0.37 

0.71 

6.7 

4.8 

5.7 

0.91 

0.20 

0.55 

6.4 

1.02 

0.04 

0.53 

4.9 

0.97 

0.65 

0.81 

7.1 

6.7 

1.27 

0.36 

0.82 

5.4 

4.9 

0.76 

0.67 

0.71 

5.8 

5.1 

5.4 

1.05 

0.32 

0.68 

MEM 

2.2 

2.9 

0.64 

0.41 

0.53 

Average 

6.1 

KOI 

5.6 

inn 

250 


550 

SAMPLE  NUMBER 
(b) 


Fig.  9  —  Channel  2  of  a  long-wave  scan  across  the  moon.  The  artifact  appearing  near 
sample  #330  in  the  uncorrected  data  (part  a)  is  largely  eliminated  in  the  corrected  data 
(part  b).  Data  obtained  at  Montauk  Point,  LI,  at  8:48  PM  on  14  August  1983. 
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FOCAL  PLANE  EQUIVALENT  CIRCUIT 

The  starting  point  for  our  analysis  is  the  preamp/ bias  circuit  diagram  and  focal  plane  equival 
circuit.  Fig.  10,  developed  with  the  aid  of  Raytheon  engineers. 


FOCAL  PLANE  GROUND 
(GOLD  COMMON) 


Fit  10  —  Bin/ input  stage  schematic  and  detector  array  equivalent  circuit 
for  long-wave  photoconductive  HgCdTe  array.  Adapted  from  correapon- 
dence  between  J.  Fattei  and  A.  Krutchkoff  (Raytheon  Miailr  Syatema). 


Note  that  the  photoconductor  Thivenin  equivalent  voltages  £,(/)  are  the  "idea T  volt 
waveforms  that  we  wish  to  recover  by  data  reduction  (Fig.  10).  The  measured  waveform  Fj,(r)  is  f 
portional  to  the  currant  /.(/)  flowing  through  the  photoconductor.  The  circuit  ground  and  focal  pi 
ground  (i.e.,  gold  common)  have  effectively  zero  resistance.  As  shown  in  Fig.  10,  the  focal  plane  a 
mon  and  circuit  ground  are  joined  by  two  resistive  ground  connections  that  an  electrically  in  para 
The  detector  equivalent  resistances  R2,  Hj,  ...Rig,  are  assumed  to  be  unknown  and  genet 
different  from  one  another.  Our  model  is  general  enough  to  accommodate  nonuniformity  in  dete 
responsivity,  accounted  for  as  an  unknown  gain  parameter  K„  multiplying  each  of  the  Thivenin  vol 
sources  in  Fig.  10. 

Figure  10  can  be  simplified  somewhat  by  replacing  the  bias/input  stage  schematics  by  equivi 
circuits.  We  note  that  the  voltage  drop  from  transistor  emitter  to  base  is  a  small  constant  value, 
points  A  and  B  in  Fig.  10  are  effectively  shorted.  Moreover,  the  large  value  of  capacitor  C(  assures 
Ci,  also,  can  be  considered  effectively  a  short  circuit  The  resultant  focal  plane/ input  stage  equivi 
circuit  is  given  in  Fig.  11.  In  this  figure  we  consider  the  currents  /.(r)  as  the  measurable  quantifier 
voltage  sources  £„(/)  as  the  quantities  we  wish  to  determine.  The  31  constants  Rf,  Jt„  am 
Or  —  2.  3,  ...  16)  are  ail  indeterminate  at  this  stage  of  analysis. 


It  follows  from  Fig.  1 1  that 


K,E,  -  RmIm  +  Rtl,  a  -  2,  3,  ...  16 


where 


'-S'- 
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Pig.  11  —  Long-wave  focal  plane/input  stage  equivalent  circuit,  corresponding  to  Fig.  10 


The  quantities  E„,  /„,  and  /  in  Eqs.  (2)  and  (3)  are  time  waveforms,  i.e.,  E„(t),  /„( t),  and  /(/).  Other 
quantities,  i.e.,  K„,  R„,  and  Rg  are  time-independent.  We  see  directly  from  Eq.  (2)  that  if  Rg  were 
zero,  the  desired  waveforms  E„(t)  would  be  equal  to  the  measured  waveforms  /„(t)  to  within  a  readily 
determined  calibration  gain.  Nonzero  Rg,  however,  causes  coupling  among  the  waveforms.  Equation 
(2)  may  be  thought  of  as  an  error  correction  equation:  if  the  constants  K„,  /?„,  and  Rg  are  known,  Eq. 
(2)  can  be  used  to  convert  the  measured  uncorrected  waveforms  /„  into  the  desired  waveforms  £„. 


Equation  (2)  is  now  written  as 
where  we  define 


AnE„  "  4  +  5,/, 
Am  m  KjRn  . 


(4) 

(5) 


and 


B„  =  Rg/R„  .  (6) 

Summing  Eq.  (4)  for  values  of  n  from  2  to  16,  using  Eq.  (3)  to  simplify  the  result,  and  gathering  terms 


in  /  we  obtain 


-rw  5 

Substituting  Eq.  (7)  into  Eq.  (4)  we  obtain 


16 


-i 


16 


AjEj. 


(7) 


4  “  A*E*  -  «„ 


16 


£  ajEj  > 


(8) 


where  we  define 


16 


=  ^|i  +  z  BA  ■ 

m-2 


-1 


(9) 


We  think  of  Eq.  (8)  as  describing  how  the  desired  waveforms,  E„(t),  are  damaged  in  creating  the 
measured  waveforms  /.(f). 


From  Eqs.  (3)  and  (4), 


E„  -  A~niH  +  B„ 


16 


(10) 


11 


■**>**»'.*.; 
*  «r,  »  ,  *  «  ' 


V  g=s  ^ g  y* . » A.  ■  .  •  A  A  .ILWL1  HW-  IU1.  1 J  J»  J  J  J,  J  .1  IPJ . .  f.  ~  J.  j.  .■  ■  j. 
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We  can  show  from  Eq.  (9)  that 


16 


(11) 


Equation  (10),  like  Eq.  (2),  may  be  thought  of  as  an  error  correction  equation.  Our  task  is  now  to 
obtain  values  for  the  30  constants  (A„,  B„)  appearing  in  Eq.  (10). 

BAR-TARGET  EXCITATION 

We  now  assume  that  the  sensor  is  scanned  across  a  vertically  oriented  bar  target,  so  that 

\EtU),  »  “  2,  4,  ...  16, 

E*U)  "  (£0(f).  n-3.  5.  ...  15.  (12) 

The  forms  of  £,(r)  and  £0(r)  are  schematized  in  Fig.  12.  The  odd-numbered  waveforms  are  delayed 
two  dwell-times  relative  to  the  even-numbered  waveforms  due  to  the  gold  common  being  equal  in 
width  to  the  detectors  (Fig.  1). 


Fig.  12  —  Idealized  even-channel  waveforms,  E,U),  and 
odd-channel  waveforms,  E0(t),  for  bar  target  excitation. 
Waveforms  E0(t)  are  delayed  two  dwell-times  relative  to 
E,(t),  i.e.,  f|  —  2t4,  due  to  the  separation  between  the 
even  and  odd  detector  columns  shown  in  Fig.  1.  Step- 
response  rise-time,  r2,  is  assumed  also  to  be  two  dwell- 
times,  including  the  effects  of  optical  blur,  detector  size, 
and  the  electrical  response  characteristics  of  preamps,  and 
digital  recorder. 


TIME 

From  Eqs.  (8)  and  (12) 


h~ 


where  we  define 


w*.*  ».  •  . 


k  -9.  -. 


.  •*  - 


v»-i 

f  ■  ■■  ■  ■■  t 

S.-V.N' 

v'  \* 


H„(Ee  -  8„£0),  n  even 

(13a) 

•  •  ;\r- 

Hh(Eq  -  8 „£*),  n  odd, 

(13b) 

PK 

'#$*:$■ 

(A„  -  «„Q),  n  even 

(14a) 

\y.;-\y 

f .  - .  -. 

CvV/S* 

,  -  €nC0),  n  odd, 

(14b) 

r..* 

i’-r  ■+ 
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and 


a.  s 


and 


'  h  even 

(15a) 

[  n  odd. 

(15b) 

re  defined  as 

c*  -  2^  Ay, 

(16) 

7 

-0-  I  A  2j+l  ■ 

(17) 

For  even-numbered  channels  the  form  of  /a(r)  follows  from  Eq.  (13a),  depicted  graphically  as 
Figs.  13a  and  13b.  Also  shown  as  Fig.  13c  is  a  more  realistic  waveform  adapted  from  Fig.  3.  The 
quantities  x,  and  shown  on  Fig.  13c  are  attributes  that  may  be  measured  directly  from  the  available 
waveforms  lnU). 


Fif.  13  —  (•)  and  (b)  The  form  of  /,(/)  it  deduced 
from  Eq.  (13a)  and  Fif.  12,  for  even-numbered  chan¬ 
nels.  (c)  A  more  realistic  waveform  adapted  from  Fig. 
3.  Equations  for  the  correction  constants  are  obtained 
by  equating  corresponding  quantities  in  (b)  and  (c). 
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Setting  corresponding  quantities  equal  to  one  another  from  Figs.  (13b)  and  (13c)  we  obtain 

.x,  -  (A  -  a)H„ 


y„  -  Mi  -  «)*„• 

Recasting  Eq.  (18), 

Hh  -  xj(b  -  a), 

and  dividing  Eq.  (19)  by  Eq.  (18), 

-  yJxH. 

Eliminating  H„  from  Eqs.  (14a)  and  (l$a)  we  obtain 

A„  ■  ca(C(  +  C(/4»)  ■  ,«)l«  '(Oi  £®*)> 
where  the  final  equality  is  obtained  by  substituting  Eq.  (21)  for  8*. 

Eliminating  Hh  from  Eqs.  (15a)  and  (20), 

Hm  —  xj(b  -  a)  -  «aC</8, |. 
Substituting  Eq.  (21)  for  8*  into  Eq.  (23)  we  can  write 

*»  “ 

where  f*  is  independent  of  n.  From  Eqs.  (22)  and  (24) 

A„  -  r,(C,y,  +  Co*). 

Summing  Eq.  (25)  over  even  values  of  n  and  invoking  Eq.  (16)  we  obtain 

c,-  r,(c,n  +  Co*). 


where  we  define 


«•> 

Eliminating  Tt  from  Eqs.  (24)  to  (26)  we  obtain 

— -  ,2,) 

and 

<„-y,C'(C,Yt  +  CQX,)-\  n  even.  (30) 

For  odd-numbered  channels  the  form  of  /„(/)  follows  from  Eq.  (13b),  depicted  graphically  as 
Figs.  14a  and  14b. 

We  can  show  analogous  to  Eqs.  (29)  and  (30)  that  for  odd-numbered  channels 

.  ~  (  Coy*  +  Ctx„  |  nn 


Ak  “  Co 


CqYq  +  c,x0 


,  h  odd. 


<7AV  •  >  "VW?  ?.  >V:; 
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Our  correction  formula  expressed  as  Eq.  (10)  has  30  indeterminate  constants,  15  values  each  of 
A„  and  B„.  However  we  see  that  if  the  A„  can  be  determined,  the  values  of  B„  can  then  be  calculated 
from  Eqs.  (35),  (16),  and  (17).  Thus,  our  30  unknowns  are  presently  reduced  to  IS.  (The  constants 
Yt,  X„  Y0,  and  X0  appearing  in  Eq.  (35)  are  calculated  from  Eqs.  (27),  (28),  (33),  and  (34)  in  terms 
of  the  x„,  y„  values  directly  measured  from  the  bar  target  waveforms,  e.g..  Figs.  3  and  4.) 

ITERATIVE  SOLUTION 

The  constants  A„  appearing  in  the  data  correction  equation,  Eq.  (10),  are  obtained  by  iteratively 
solving  Eqs.  (29)  and  (31)  according  to  the  following  prescription: 


A 


where,  from  Eqs.  (16)  and  (17), 


1  +  C^x„  1 

|  cf0  rt  +  c0(nxe )’  neven 

(36) 

1  Ct‘>yM  +  CtU)x,  ) 

1  C0(,)  Yo  +  CtU)X0  (’  ',0dd’ 

(37) 

c"  -  i;  Atf. 

(38) 

Co(,)  -  TAjjl,. 

(39) 

The  superscript  on  quantities  A„,  C, ,  and  C0  in  Eqs.  (36)  to  (39)  indicates  the  order  of  iteration. 

We  start  the  iterative  solution  by  assuming  that  to  zero-order  the  gains  K„  and  resistances  R„  are 
the  same  for  all  detectors  (cf.  Fig.  10).  From  Eq.  (5), 

Aim-l.  (40) 

Starting  with  Ai0)  equal  to  a  constant  other  than  unity  results  in  a  proportional  scaling  of  the  iterative 
solution  for  A„.  Since  the  measuiables  V„(t)  in  Fig.  10  are  related  to  the  /„(/)  by  an  electrical  gain,  it 
follows  from  Eq.  (8)  that  all  values  of  A„  can  be  scaled  by  an  arbitrary  multiplicative  constant  chosen  to 
suit  our  convenience. 

From  Eqs.  (38)  to  (40), 

C®>  -  8  (41) 

Ci0)  -  7.  (42) 

The  next  order  of  approximation  beyond  Eq.  (40)  follows  from  Eqs.  (36),  (37),  (41),  and  (42), 

*  (»8£ +??.)•  "even 

A<“-  («) 


Ty*  +  *x, 
7Y0  +  8X0  ’ 


n  odd. 


The  parameters  x„  and  y„  in  Eq.  (43)  are  directly  measurable  from  the  bar  target  data,  e.g.,  Figs.  3  and 

The  numerical  values  A,*0  obtained  from  Eq.  (43)  are  substituted  into  Eqs.  (38)  and  (39)  to 
obtain  C,(1)  and  Cj°,  which  in  turn  are  substituted  into  Eqs.  (36)  and  (37)  to  obtain  A„<v.  Conver¬ 
gence  is  achieved  to  eight  places  after  four  iterations. 
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The  fully  converged  solution  for  A„  is  used  with  Eqs.  (16),  (17),  and  (35)  to  obtain  the  remain¬ 
ing  correction  constants,  B„. 

The  waveform  parameters  x„,  y„  obtained  from  the  Fig.  2  waveforms  and  the  corresponding  itera¬ 
tive  solutions  for  A~'  and  B„  are  given  in  Table  2. 

CONCLUSIONS 

As  shown  in  the  numerical  examples  discussed  in  our  Summary  of  Results,  illustrated  in  Figs.  2 
to  9,  excellent  correction  has  been  achieved  for  data  defects  caused  by  resistive  interchannel  coupling. 
The  correction  constants  derived  from  a  "diagnostic  frame,"  Fig.  2,  provide  accurate  correction  to 
scenes  containing  different  varieties  of  structure  and  obtained  months  after  the  diagnostic  frame.  Thus, 
the  correction  constants  are  time-  and  data-independent.  Nonetheless,  we  intend  to  obtain  diagnostic 
frames  both  immediately  before  and  immediately  after  each  field  trial  with  the  Raytheon  sensor  to 
assure  a  completely  updated  set  of  long-wave  resistive  coupling  correction  constants  for  use  in  data 
reduction. 

The  instantaneous,  "memoryless,"  nature  of  the  resistive  coupling  defect  allows  us  to  perform 
correction  without  smoothing  the  data,  i.e.,  without  loss  of  spatial  resolution.  The  form  of  the  correc¬ 
tion  equation,  Eq.  (10),  together  with  the  fact  that  B„  -  O(10"2),  implies  that  the  relative  increase  in 
rms  sensor  noise  introduced  by  the  correction  process  is  a  small  fraction  of  1%,  and  hence  is  negligible. 
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Appendix 

DERIVATION  OF  EQUATION  (35) 


From  Eqs.  (30)  and  (27) 


£ 


«2.-c(n(c,n  +  w‘. 


From  Eqs.  (32)  and  (33) 


jr  <2n+i  -  c0ro(c0r0  +  Qat0)-«. 

Adding  Eqs.  (Al)  and  (A2),  and  subtracting  the  resultant  from  unity,  we  can  show  that 


v  CtC0(XtX0-  YeY0) 

*"  ”  (C0  Jo  +  C'X o)  (C,  n  +  C0*,)  ' 


From  Eqs.  (A3),  (30),  and  (32), 


it 


£«.|  - 

m—2 


F0  +  (C'/CQ)X0 1 

,,even 

n  +  (c0/ce)xt  i 


IV  Jo 


,  n  odd. 


Equation  (A4)  is  given  as  Eq.  (35)  in  the  main  text. 
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